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Abstract 


The  tunneling  effects  in  the  molecular  inversion  of  aziridine  are  investigated.  A 
full-dimensional  potential  was  constructed  based  on  ab  initio  results  and  used  to  study  the  effect 
of  vibrational  excitations  on  tunneling.  Using  the  semiclassical  approach  that  incorporates 
tunneling  into  classical  trajectory  calculations,  it  is  found  that  excitations  of  all  the  modes  except 
the  reaction  coordinate  have  negligible  effects  on  tunneling.  This  allows  accurate  thermal-rate 
calculations  carried  out  by  using  a  one-dimensional  (1-D)  model.  (This  is  in  accord  with  the 
conclusions  of  the  work  by  N.  Rom,  V.  Ryaboy,  and  N.  Moiseyev  [Chemical  Physics  Letters, 
vol.  204,  p.  175, 1993]  and  by  Z.  Smedarchina,  W.  Siebrand,  and  M.  Z.  Zgierski  [Journal  of 
Chemical  Physics,  vol.  102, p.  7024, 1995].)  A  1-D  Weutzel-Kramers-Brillouin(WKB)  formula 
was  employed  to  obtain  the  level  splitting  and  was  found  to  be  very  accurate  based  on 
comparisons  with  the  quantum  results.  Furthermore,  the  calculated  thermal  rates  are  in  good 
agreement  with  the  experimentally  measured  values.  The  results  provide  insight  into  the  reaction 
mechanism  and  explanations  for  the  experimental  findings  of  non-Rice-Ramsperger- 
Kassel-Marcus  (RRKM)  behavior  and  incoherent  tunneling. 
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1.  Introduction 


The  intramolecular  conversion  in  aziridine  (see  Figure  1)  has  been  extensively  studied  both 
experimentally  and  theoretically  [1-5].  Of  particular  interest  is  the  finding  by  Borchardt  and 
Bauer  [1]:  the  observed  inversion  rate  constants  do  not  depend  on  pressure  over  the  range  of 
5-1,100  torr  for  temperatures  from  298-388  K,  while  standard  Rice-Ramsperger-Kassel-Marcus 
(RRKM)  calculations  indicate  that  the  system  should  be  in  the  bimolecular  regime  for  pressures 
below  100  torr.  This  non-RRKM  behavior  was  attributed  to  the  slow  internal  vibrational 
relaxation  (IVR)  in  this  molecule — the  rate  of  reaction  is  limited  by  the  rate  of  IVR  [1]. 
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Figure  1.  Schematic  Illustration  of  Molecular  Inversion  in  Aziridine. 

The  purpose  of  this  paper  is  to  provide  a  comprehensive  study  of  the  inversion  dynamics  of 
aziridine  to  gain  insight  about  the  reaction  mechanism.  Of  main  interest  is  in  the  tunneling 
effects  in  this  reaction.  Since  the  migrating  atom  is  hydrogen  and  the  experimental  temperatures 
are  near  300  K,  appreciable  tunneling  effects  are  expected.  The  method  employed  here  is  the 
multidimensional  semiclassical  approach  that  has  been  used  for  unimolecular  and  bimolecular 
reactions  [6-13],  which  incorporates  tunneling  within  classical  trajectory  simulations  and,  thus, 
allows  for  explicit  treatment  of  the  full-dimensional  dynamics. 

There  are  already  some  theoretical  studies  on  tunneling  in  aziridine  [4,  5]  in  which  a 
one-dimensional  (1-D)  potential  along  the  reaction  coordinate  [4]  or  an  adiabatic  1-D  potential 
that  includes  the  effects  of  high-frequency  modes  [5]  is  used.  Rom,  Ryaboy,  and  Moiseyev  [4] 
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treated  the  effect  of  other  degrees  of  freedom  by  introducing  an  absorbing  potential.  They  also 
performed  calculations  using  a  two-dimensional  (2-D)  potential,  which  includes  the  N-H  stretch 
mode,  and  obtained  good  agreement  with  the  results  from  the  1-D  calculations.  Both  studies 
[4,5]  yielded  good  agreement  with  the  experimental  thermal  rates,  suggesting  that  a  1-D 
treatment  is  adequate.  The  multidimensional  semiclassical  approach  enables  direct  verification 
of  this  point  by  explicitly  investigating  the  dynamic  motions  of  all  the  atoms  and  their  effects  on 
tunneling. 

The  physical  quantity  of  interest  is  the  thermal-rate  constant,  for  which  experimental  results 
for  both  the  normal  and  deuterated  species  are  available  [2].  In  principle,  it  can  be  calculated  by 
the  standard  expression 


v  k°exp(-En/kBT) 

Inexp(-En/kBT) 


(1) 


where  En  is  the  total  energy  of  the  molecule  in  a  specified  state  with  quantum  numbers 
n  =  {m,  n2. . .  >,  and  kn  is  the  reaction  probability  for  that  state.  However,  the  number  of  terms  in 
the  sum  can  be  quite  large  as  the  energy  increases,  making  the  computation  prohibitively 
expensive.  With  the  use  of  the  multidimensional  semiclassical  approach,  it  is  explicitly 
demonstrated,  by  selectively  exciting  each  mode,  that  there  is  only  one  mode  that  plays  a 
dominant  role  in  the  tunneling.  Thus,  the  thermal-rate  calculation  reduces  to  a  much  simpler  1-D 
problem.  Since  the  reaction  probability  for  a  given  energy  would  be  independent  of  die 
excitation  site  if  the  system  were  statistical,  these  mode-selective  results  also  provide  insight  into 
the  observed  non-RRKM  behavior. 

Besides  the  apparent  non-RRKM  behavior,  another  interesting  aspect  of  the  system  is  the 
incoherent  tunneling,  (i.e.,  the  reaction  is  characterized  by  rate  processes).  For  a  symmetric 
double-well  potential,  it  is  usually  expected  that  the  molecular  motions  are  coherent  oscillations 
between  the  wells  rather  than  rate  processes  and  that  the  oscillating  frequency  is  proportional  to 
the  level  splitting.  A  rate  process  occurs  in  a  double-well  system  only  when  there  is  some 
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damping  mechanism  that  broadens  the  energy  levels.  If  this  broadening  exceeds  the  level 
splitting,  then  tunneling  becomes  irreversible.  In  the  present  case,  since  the  measured  rates  [1] 
show  no  pressure  dependence,  the  intermolecular  collisions  should  not  be  the  major  cause  of  the 
broadening.  The  other  possible  cause,  as  suggested  by  Smedarchina,  Wiebrand,  and  Zgierski  [5], 
is  the  intramolecular  couplings;  the  density  of  states  and  the  anharmonicity  increase  rapidly  as 
the  energy  increases,  which  may  lead  to  greater  broadening  than  splitting  for  high  energy  states. 
The  results  here  are  in  agreement  with  those  of  Rom,  Ryaboy,  and  Moiseyev  [4]:  the  level 
splittings  of  the  low-energy  states  are  too  small  to  be  observed,  and  the  states  contributing  to  the 
measured  thermal  rates  are  those  near  the  barrier  top.  Thus,  rate  constants  rather  than  level 
splittings  were  observed  in  the  experiments  [1,2]. 

This  paper  is  organized  as  follows.  The  potential  energy  surface  is  described  in  section  2,  the 
semiclassical  method  for  computing  thermal  rates  in  section  3,  the  results  and  discussions  in 
section  4,  and  the  conclusions  in  section  5. 

2.  Potential  Energy  Surface 

The  potential  energy  surface  used  in  this  study  is  based  on  a  valence  force  field  with  the 
normal-mode  frequencies  fitted  to  ab  initio  results  calculated  using  the  Gaussian  94  codes  [14]. 
Both  the  equilibrium  and  transition  states  are  calculated  at  the  MP2  [15]  level  with  the  frozen 
core  approximation  using  the  6-311G/(2df,  2p)  [16-19]  basis  set.  The  functional  form  of  the 
total  potential  is 


v = S,  (♦)’ VA  +  [l  -  S,  (♦)]' VB  + Vfo).  (2) 

Here,  <|>  is  the  reaction  coordinate  (the  NH  wag  angle),  V(<}>)  is  a  double-well  potential  along  <|>, 
Va  and  Vb  consist  of  the  terms  corresponding  to  the  two  isomers,  and  Si(<]>)  is  a  switching 
function  that  smoothly  connects  Va  and  Vb,  that  is. 
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Si  (<!>)=  l[l  +  tanh(a<|>)]. 


(3) 


where  a  is  chosen  to  be  3.0  rad’1. 

The  double-well  potential  V(<j>)  is  of  the  form 

V((())=C0+C2(l)2+C4(t>4,  (4) 

where  Co  =  19.2kcal/mol  and  is  the  barrier  height,  and  C2  and  C4  are  27.1  and  9.57 
rad_4kcal/mol,  respectively,  determined  by  requiring  V(<J>)  and  its  first  derivative  at  the 
equilibrium  NH  wag  angle  (68.2°)  to  be  zero. 

Each  of  the  potentials,  Va  and  VB,  is  a  sum  of  Morse  functions  for  bond  stretches,  harmonic 
functions  for  the  angle-bending  terms,  and  three-term  cosine  expansions  for  the  dihedral 
interactions,  that  is, 

^A,  B  —  ^stretch  +  Vbend +VdaKdral ,  (5) 

where 

Vstretch  =SDi  {-exp  [- P (rs  -r° )]}2 ,  (6) 

l 

v^=Iik  ,(e,-0»f,  (7) 

1 

and 

^dihedral  =?k  +  ai2  cos  (2t)+  ai4  cos  (4t)].  (8) 

i 

The  values  of  the  parameters  in  equations  (6) — (8)  are  given  in  Table  1 .  The  equilibrium  CNH 
angles  used  in  the  bending  potential  are  also  functions  of  the  reaction  coordinate 
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Table  1.  Potential  Energy  Surface  Parameters 


Morse  Potential  V(r) 

Bond 

r° 

D 

p 

0 

(A) 

(kcal/mol) 

(A'1) 

N-H 

1.014 

99.77 

2.135 

C-C 

1.478 

54.73 

1.966 

C-N 

1.473 

45.89 

2.678 

C-H' 

1.079 

110.61 

1.826 

C-H" 

1.080 

110.61 

1.825 

Harmonic  Bending  Interactions  V(0) 

Angle 

0 

0 

k 

(°) 

(kcal/mol  rad2) 

HNC 

108.75 

74.32 

NCC 

59.90 

113.45 

CNC 

60.21 

183.18 

NCH" 

118.32 

77.78 

NCH' 

114.33 

99.61 

CCH" 

117.77 

94.23 

CCH' 

119.53 

100.89 

Torsional  Interactions,  V(Tj) 

Dihedral  Angle 

T° 

ao 

a2 

a4 

(°) 

(kcal/mol) 

(kcal/mol) 

(kcal/mol) 

H'CCH" 

149.04 

1.738 

-2.267 

1.204 

9(<>)=9„-(e„  eJexpl-bV),  (9) 

where  Qeq  =  108.8°,  0ts  =  149.8°,  and  b  is  chosen  to  be  5.0  rad'4. 


The  comparison  of  the  scaled  (with  a  scaling  factor  0.95)  ab  initio  and  die  fitted  frequencies 
is  given  in  Table  2.  It  is  noted  that  an  accurate  fitting  is  not  crucial  here  since  the 
full-dimensional  force  field  is  only  used  in  a  qualitative  mode-selective  study. 
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Table  2.  Normal-Mode  Frequencies 


Mode 

Ab  Initio 
(cm'1) 

Analytical  PES 
(cm"1) 

7 

745 

699 

8 

821 

788 

9  1 

842 

814 

10 

876 

853 

11 

967 

957 

12 

1,065 

1,067 

13 

1,068 

1,073 

14 

1,108 

1,136 

15 

1,188  ~1 

1,194 

16 

1,203 

1,221 

17 

1,242 

1,738 

18 

1,433 

1,441 

19 

1,459 

1,396 

20 

3,018 

3,011 

21 

3,024 

3,017 

22 

3,108 

3,113 

23 

3,119 

3,124 

24 

3,380 

3,380 

3.  The  Semiclassical  Method 


The  semiclassical  method  for  computing  tunneling  splitting  has  been  described  previously 
[9, 13].  The  new  feature  here  is  that  a  different  Weutzel-Kramers-Brillouin  (WKB)  expression 
for  level  splitting  in  a  1-D  potential  is  used.  This  more  complicated  form  is  valid  for  the  entire 
range  of  energies  and  reduces  to  the  one  used  previously  for  low  energies.  It  is  of  the  form  [20] 


AE  =  4/iv 


Ke-25-'*  +ljn -1 

l|(e-2s-,’,+l)'!+r 


00) 


where  v  is  the  frequency  of  the  oscillator  at  energy  E  and  Sc  is  the  classical  action  integral 
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(11) 


Sc  =  fs>  V2m[v(s)-E]ds. 

Using  this  WKB  expression  for  a  1-D  potential,  the  splitting  for  a  multidimensional  system 
can  be  approximated  by  [9, 13] 


AE  =  2»l<A(t), 

with  A  (t)  accumulated  by  propagating  a  classical  trajectory 


(12) 


A(t)=2Eh(“°JP^ir^-  (13) 

Here,  h(t  -  tn)  is  the  usual  step  function  and  tn  are  the  times  that  a  trajectory  is  at  the  turning 
points.  The  brackets  in  equation  (12)  imply  an  average  over  an  ensemble  of  trajectories  with  the 
same  quantum  number  but  different  vibrational  phases. 

The  key  issue  is  the  evaluation  of  the  action  integral  Sc.  It  was  found  in  previous  studies 
[4, 5]  that  the  only  mode  strongly  coupled  to  the  reaction  coordinate  in  the  barrier  region  is  the 
N-H  stretch — the  effect  of  which  results  in  a  reduction  of  the  effective  barrier  height  by  about 
1  kcal/mol.  Considering  the  uncertainty  in  the  ab  initio  calculation  of  the  barrier  height  and 
width,  it  is  reasonable  to  neglect  the  intramolecular  couplings  and  just  use  the  potential  along  the 
reaction  coordinate  V(<]>)  [given  by  equation  (4)]  to  compute  the  action  integral.  It  is  noted  that 
the  turning  points  are  still  determined  in  the  full-dimensional  space  by  propagating  trajectories. 
The  1-D  potential  is  used  only  in  evaluating  the  action  integral. 

The  mode-selective  results  described  in  the  next  section  show  that  excitation  of  the  reaction 
mode  gives  the  dominant  contribution  to  tunneling.  Thus,  in  computing  the  thermal  rate,  it  is 
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assumed  that  all  the  modes  except  the  tunneling  coordinate  remain  in  their  ground  states  and  the 
thermal  rate  is  approximated  by 


,/tn  EnkDexp(-ea/kBT) 

[  )=  Enexp(-£n/kBT)  ’ 


(14) 


with  £n  being  the  energy  level  in  the  1-D  potential  V(<t>).  The  reaction  rate  kn  is  given  by 


kn=vn 


(15) 


for  En  >  Vo,  and 


kn=rn//*=-fe£ 


4 aV 


(16) 


for  £„  <  Vo.  Here,  vn  is  the  oscillating  frequency  at  energy  £n,  V0  is  the  barrier  height  of  V(<j>),  Tn 
is  the  level  width,  and  AEn  is  the  level  splitting.  The  last  equation  can  be  seen  from  the  WKB 
expressions  for  the  level  width  and  splitting  at  low  energies 

T = ^v  exp  (-  2SC  !h)  (17) 


and 


AE=  2fivexp(-Sc/ft)  .  (IB) 

4.  Results  and  Discussion 

Mode-selective  effects  in  aziridine  were  investigated  first  using  the  full-dimensional 
potential.  The  calculations  were  carried  out  for  initial  conditions  with  5  kcal/mol  of  energy 
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assigned  to  each  normal  mode  above  the  zero-point  level.  The  results  show  that  excitation  of  the 
reaction  coordinate  mode  has  the  dominant  effect  on  tunneling,  with  the  splitting  being  at  least 
three  orders  of  magnitude  larger  than  that  for  excitation  of  any  other  mode.  This  provides  some 
insight  into  the  reaction  mechanism  and  the  observed  non-RRKM  behavior.  The  mode-selective 
results  indicate  that  the  intramolecular  couplings  may  be  weak  enough  that  the  rate  of  IVR  is  the 
slower  and,  hence,  rate-limiting  step.  Thus,  the  observed  reaction  rates  are  independent  of 
pressure,  and,  since  excitations  of  many  of  the  modes  have  no  contributions  to  the  state  counting, 
the  observed  rates  at  high  pressures  deviate  from  those  of  RRKM  predictions  (see  Figure  6  of 
Borchardt  and  Bauer  [1]). 

Based  on  the  mode-selective  results  and  the  fact  that  the  couplings  between  the  reaction 
coordinate  and  the  rest  of  the  modes  are  not  strong  enough  to  significantly  affect  the  effective 
barrier,  it  is  reasonable  to  expect,  as  suggested  previously  [4, 5],  that  a  1-D  treatment  is  adequate. 
The  thermal  rate  is  thus  computed  using  equation  (14)  for  the  1-D  potential  along  the  reaction 
coordinate,  V(<|>).  From  equations  (15)  and  (16),  it  is  clear  that  one  needs  to  obtain  the  oscillating 
frequency  vn  and  the  level-splitting  AEn  for  all  the  states  contributing  to  the  thermal  rate.  The 
frequency  vn  was  computed  by  first  using  adiabatic  switching  to  confine  the  trajectories  to  the 
semiclassical  eigen  state  and  then  propagating  them  one  cycle  to  determine  the  period.  The 
level-splitting  was  obtained  from  equation  (10)  by  using  the  adiabatically  switched  trajectories  to 
compute  the  action. 

It  is  important  to  examine  the  accuracy  of  the  semiclassical  WKB  formula  equation  (10). 
Table  3  shows  comparisons  of  the  semiclassical  and  exact  quantum-mechanical  level  splittings 
for  both  the  normal  and  deuterated  species.  The  semiclassical  splittings  were  computed  by  using 
equation  (10),  and  the  quantum  ones  were  obtained  by  standard  diagonalization  of  the 
Hamiltonian.  The  semiclassical  and  quantum  results  agree  to  within  a  few  percent, 
demonstrating  the  accuracy  of  the  WKB  expression.  (The  quantum  splittings  for  the  first  three 
deuterated  states  are  missing  because  the  accuracy  of  the  calculation  is  not  high  enough.) 
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Table  3.  Quantum  Mechanical  and  Semiclassical  Tunneling  Splittings. 


n 

Hydro) 

?en 

Deuterium 

Quantum  Mechanical 
(kcal/mol) 

Semiclassical 

(kcal/mol) 

Quantum  Mechanical 
(kcal/mol) 

Semiclassical 

(kcal/mol) 

0.458  x  10' 12 

— 

0.682  x  10' 18 

1 

SimS 

wsebesM 

— 

0.354  x  10" 15 

2 

0.248  x  10"7 

— 

0.828  x  10' 13 

3 

0.120  x  10"10 

0.118  x  10' 10 

4 

IKHmiSHi 

TTTHM 

0.116  x  10'8 

0.115  x  10"8 

5 

6 

0.924  X  10'1 

0.953  x  10'1 

0.420  x  10"5 

HEE333HI 

7 

— 

— 

■buheshh 

8 

— 

— 

1 

|  9 

— 

— 

0.707  X  10'1 

0.731  x  10' 1  | 

Figure  2  provides  a  comparison  of  the  calculated  and  experimental  thermal-rate  constants  for 
both  the  normal  and  deuterated  species.  The  barrier  height  was  adjusted  to  19.2  kcal/mol  to  fit 
the  experimental  result  for  hydrogen.  This  leads  to  good  agreement  for  deuterium.  The  adjusted 
barrier  height  is  close  to  the  ab  initio  value  of  18.2  kcal/mol,  indicating  that  the  intramolecular 
couplings  do  not  significantly  contribute  to  the  effective  barrier. 


To  illustrate  the  contributions  from  both  classical  and  tunneling  regimes,  the  contributions  to 
the  thermal  rate  of  each  individual  state,  kn  exp  (-£n/kBT),  are  plotted  in  Figure  3.  The  two 
curves  correspond  to  T  =  300  and  360  K,  which  are  near  the  low  and  high  ends  of  the 
experimental  data  shown  in  Figure  2.  The  barrier  height  is  slightly  below  87.  It  is  clear  that,  for 
the  temperature  range  of  the  experiments,  the  major  contributions  to  the  rate  are  from  the  three 
states  near  the  barrier:  one  below  (n  =  6)  and  two  above  (n  =  7,  8).  Thus,  tunneling  contributes 
about  1/3  of  the  total  rate.  Table  4  provides  the  quantitative  results  on  the  tunneling  correction 
factor  as  a  function  of  temperature  defined  as 


^  ^  (T ) total  /  ^  (T  Classical 


(19) 
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Figure  2.  The  Thermal  Rates  as  a  Function  of  Temperature  for  Both  the  Normal 

and  Deuterated  Species.  The  Points  and  Lines  Represent  Experimental  and 
Calculated  Results,  Respectively.  The  Experimental  Data  Are  Taken 
From  [2], 


Figure  3.  The  Contribution  to  the  Thermal  Rate  From  Each  Individual  State, 
k„exp  (-ej/kfiT).  The  Barrier  Is  Slightly  Below  £7. 

The  values  for  both  the  normal  and  deuderated  species  are  given.  The  results  obtained  by  using 
the  Wigner  formula  [21] 


KWigner 


1 

barrier 

24 

kT 

(20) 
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Table  4.  Tunneling  Correction  Factors  Calculated  by  Using  the  1-D  Potential  V  (<J>). 
The  Results  Are  From  the  Exact  Calculation  Employing  Equation  (19) 
and  by  Using  the  Wigner  Formula  Equation  (20) 


■ 

Hydrogen 

Deuterium 

K 

^Wigner 

K 

^Wigner 

260 

1.88 

1.79 

1.32 

1.40 

280 

1.66 

1.68 

1.26 

h  1.34 

300 

1.52 

1.59 

1.22 

1.30 

•320 

1.42 

1.52 

1.19 

1.26 

340 

1.35 

1.46 

1.17 

1.23 

360 

1.30 

1.41 

1.15 

1.21 

380 

1.25 

1.37 

1.13 

1.18 

400 

1.22 

1.33 

1.12 

1.17 

are  also  given  for  comparison.  It  is  seen  that  the  Wigner  formula  works  well,  despite  the  fact 
that  it  is  derived  for  the  single  barrier  (no  well)  case,  while  this  case  is  a  double-well  system. 


At  low  energies  where  the  energy  levels  are  widely  separated,  the  system  should  exhibit 
periodic  oscillations  between  the  wells,  with  the  oscillating  frequency  proportional  to  the 
tunneling  splitting,  which  can  be  measured  spectroscopically.  At  high  energies  where  the 
density  of  states  is  high  and  the  intramolecular  couplings  are  strong,  tunneling  becomes  a  rate 
process  when  the  level  broadening  exceeds  the  splitting.  These  calculations  provide  an 
explanation  for  experimental  observation  of  rate  processes  rather  than  level  splittings  [1, 2];  the 
splittings  of  the  low-lying  states  are  too  small  to  be  measured;  the  contribution  to  the  thermal 
rates  comes  mainly  from  the  states  near  the  barrier  top,  where  the  anharmonicity  is  large  and  the 
density  of  states  is  high  due  to  the  large  number  of  degrees  of  freedom  of  the  molecule;  hence, 
reaction  rates,  rather  than  splittings,  were  observed. 


5.  Conclusions 


The  intramolecular  conversion  in  aziridine  has  been  investigated.  Ab  initio  calculations  at 
the  MP2  level  have  been  performed,  and  the  results  have  been  used  to  construct  a  potential 
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energy  surface  that  includes  all  18  vibrational  degrees  of  freedom.  This  full-dimensional 
potential  was  used  to  study  the  mode-selective  effects  by  employing  the  multidimensional 
semiclassical  approach  [13],  and  it  was  found  that  excitations  of  modes  other  than  the  reaction 
mode  make  negligible  contributions  to  the  thermal  rate.  Thus,  the  previous  conclusion  [4,  5]  is 
verified  directly,  that  is,  the  problem  of  computing  thermal  rates  can  be  reduced  to  a  simple  1-D 
case. 

This  study  provides  some  explanations  for  the  experimental  findings  [1,  2].  The 
mode-selective  results  indicate  that  the  intramolecular  couplings  may  be  weak,  and,  thus,  the 
slow  IVR  is  the  rate-limiting  step,  leading  to  the  pressure-independent  reaction  rates.  Moreover, 
the  observed  rates  at  high  pressures  deviate  from  those  of  RRKM  calculations  [1]  because  many 
of  the  modes  make  no  significant  contributions.  The  results  also  show  that  the  level  splittings  of 
the  low  energy  states  are  too  small  to  be  observed,  and,  in  the  experimental  temperature  range, 
the  states  near  the  barrier  top  dominate  the  rate.  These  are  consistent  with  the  results  of  Rom, 
Ryaboy,  and  Moiseyev  [4].  Thus,  rate  processes  rather  than  splittings  were  observed  due  to  large 
anharmonicity  and  high  density  of  states  near  the  barrier  top. 


Intentionally  left  blank. 
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